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Abstract 

The helicity structure of the diffractive electroproduction of p mesons, e+p — ► e+p+Y, 
is studied in a previously unexplored region of large four-momentum transfer squared at 
the proton vertex, t: < t' < 3 GeV 2 , where t 1 = \t\ — \t\ m i n - The data used are 
collected with the HI detector at HERA in the kinematic domain 2.5 < Q 2 < 60 GeV 2 , 
40 < W < 120 GeV. No t dependence of the Tqq spin density matrix element is found. 
A significant t dependent helicity non-conservation from the virtual photon to the p meson 
is observed for the spin density matrix element combinations Tq + 2r\ l and Tq + 2r\± . 
These t dependences are consistently described by a perturbative QCD model based on the 
exchange of two gluons. 
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1 Introduction 



Measurements of exclusive vector meson (VM) production in ep scattering at high energy: 

e + p^e + VM + Y , (1) 

have led to considerable recent progress towards an understanding of diffraction in terms of 
QCD [jl]-@] . The reaction is induced by a real or virtual photon and Y is either a proton ("elastic" 
scattering) or a baryonic system of mass My which is much lower than the photon-proton centre 
of mass energy W ("proton dissociative" scattering). Particularly sensitive tests of QCD models 
are provided by the study of the helicity structure of the interaction and its t dependence, t being 
the square of the four-momentum transfer from the incident proton to the scattered system Y. 
The scope of this paper is to test diffractive dynamics through the extension of helicity amplitude 
extractions for exclusive p° electroproduction to larger values of \t\ than has previously been 
possible. 

Three angles are defined to characterise the electroproduction of vector mesons decaying 
into two charged particles: $ is the angle between the VM production plane (defined as the 
plane containing the virtual photon and the VM directions) and the electron scattering plane in 
the (7*1?) centre of mass system, 9* and cp are the polar and the azimuthal angles, respectively, 
of the positively charged decay particle in the VM rest frame, the quantisation axis being taken 
as the direction opposite to that of the outgoing Y system. In this paper, the distributions of the 
angles $ and 9* are analysed. 

The angular distributions give access to spin density matrix elements, which are bilinear 
combinations of the helicity amplitudes T\ VM \ , where \vm (A 7 ) is the vector meson (virtual 
photon) helicity [Q] . In the case of vector meson electroproduction by unpolarised beams and 
their subsequent decay into two pseudoscalar particles (p — > , <p — > K + K~), the $ and 9* 
distributions, integrated over the other two angles, are related to 5 of the 15 spin density matrix 
elements r\,- and r^j through the relations [|5|] 

^ oc 1 + y/2e(l + e) cos $ (rjj + 2r 5 n ) - e cos 2$ (r^ + 2rJ x ) (2) 



da 



dcos 9* 



ocl-r ° 4 +(3r ° 4 -l)cos 2 r , (3) 



where e is the polarisation parameter, i.e. the ratio of the longitudinal to transverse virtual 
photon fluxes. For this analysis e ~ 0.99 [|IJ]. 

Assuming natural parity exchange to hold (T_ Ai/m -a 7 = (— l) XvM ~ x i T Avm a 7 ), these five 
spin density matrix elements are related to the five independent complex helicity amplitudes by 
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the following relations: 



r°o 4 ocl(|T 00 | 2 + |T 01 | 2 ) 

r oo oc — Re (T 00 iji) 

r H oc i (Re (T 10 T n ) - Re (T^)) 

1 _ ^ It> 12 

r oo oc — |T 01 | 
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with N=\T 00 \ Z + \T n \' + \T 01 \* + 2 |T 10 r + \Ti-i\' . (4) 

In the case of ^-channel helicity conservation (SCHC), Xvm — A 7 , on ly the ^oo an d T n helicity 
"non-flip" amplitudes are non-zero, r^ = r n = Tq = r\ x = and the ratio i? of the cross 
sections for longitudinal to transverse photons is given by R = 1/e r^/(l — t-qq). 

In recent years, the spin density matrix elements describing process ([I]) have been measured 
for the elastic electroproduction of p and 4> mesons in the kinematic range Q 2 > 2.5 GeV 2 
and \t\ < 0.5 GeV 2 [|Tj-0], Q 2 being the negative square of the virtual photon four-momentum. 
Three main features have emerged from these measurements: 



• the dominance of the longitudinal T o over the transverse Tn helicity non-flip amplitudes; 

• the presence of a small but significant violation of SCHC, observed through the non- 
zero value of the ro matrix element, in which the dominant helicity single flip amplitude 
describes the transition from a transverse photon to a longitudinal vector meson (T i); 

• values compatible with zero for the other amplitudes describing single (T 10 ) or double 
helicity flip (T x _i). 

These features are in agreement with calculations based on perturbative QCD (pQCD) 
p. In these approaches, vector meson electroproduction is described in the proton rest frame 
as the convolution of a virtual photon fluctuation into a qq pair at a long distance from the 
target, a hard interaction mediated by the exchange of two gluons (each of them must carry 
sufficiently large transverse momentum to resolve the qq pair and the proton structure), and 
the subsequent recombination of the quark pair into a vector meson. For massless quarks, the 
helicity of the qq pair is zero, such that the helicity of the virtual photon is transferred into 
the projection of the orbital angular momentum of the qq pair onto the 7* direction. During 
the interaction, the helicity and the impact parameter of the quark pair are unchanged, but the 
orbital angular momentum can be modified through the transfer of the transverse momentum 
carried by the gluons. The helicity of the outgoing vector meson can thus be different from 
that of the incoming photon. Calculations show that such a helicity flip between the photon and 
the vector meson requires an asymmetric sharing of the photon longitudinal momentum by the 
quark and the antiquark . 

For Q 2 above a few GeV 2 and \t\%> Q 2 , the following features are expected for the (predom- 
inantly imaginary) amplitudes: 
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• a ratio constant with t for the helicity conserving amplitudes \T n \ / \T 00 \; 

• a y/\t\ dependence for the ratio of the single helicity flip to the non-flip amplitudes 

|T 01 | / |T 00 | and |T 10 | / |T 00 |; 

• a dependence linear with t for the ratio of the double flip to the non-flip amplitudes 

l^i-i I / l^ool; 

• the hierarchy: 

|T 00 | > \T n \ > \T i\ > \T W \ > . (5) 

These features are expected to hold for proton dissociative as well as for elastic scattering. 

Compared to previous results the present paper extends considerably the t range of 

the measurement of spin density matrix elements for p meson diffractive electroproduction 

e + p-^e + p + Y; p — > ir + n~ , (6) 

where the p mass range is defined by restricting the invariant mass M n7r of the decay pions to 
the interval 

0.6 < < 1.1 GeV . (7) 

Elastic and proton dissociative data are combined, and the kinematic domain of the mea- 
surement is: 

2.5 < Q 2 < 60 GeV 2 
40 < W < 120 GeV 
0< t' <3GeV 2 . (8) 

The variable t' = \t\ — \t\ min is used for the analysis, where \t\ min is the minimal value of |t| 
kinematically required for the vector meson and the system Y to acquire their effective mass 
through longitudinal momentum transfer. The t' variable, which is very well approximated 
as the square of the transverse momentum of the scattered system Y, describes the transverse 
momentum transfer to the target and is thus the relevant dynamical variable. In the elastic case 
and for moderate Q 2 , \t\ min is negligible and t' ~ \t\. 

The large t' domain covered by the present data allows for the first time a detailed study of 
the t' dependence of the helicity structure of diffractive vector meson electroproduction. 

2 Experimental procedure 

2.1 Event selection, kinematic variables and Monte Carlo simulations 

The data used for the present analysis were taken with the HI detector in 1997. The energies 
of the HERA proton and positron beamsQ were 820 GeV and 27.5 GeV, respectively. The 

1 In the following, the word electron will be used for both electrons and positrons. 
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integrated luminosity used for the analysis amounts to 6.0 pb _1 . The relevant parts of the 
detector, for which more details can be found in are the central tracking detector, the 

liquid argon (LAr) and the backward electromagnetic (SPACAL) calorimeters and the forward 
detectors, which are sensitive to energy flow close to the outgoing proton direction^ i.e. the 
proton remnant tagger (PRT) and the forward muon detector (FMD). 

Events corresponding to reaction (Q), in the kinematic range defined by relations (|]), are 
selected by requesting the reconstruction of a cluster in the SPACAL calorimeter with energy 
larger than 17 GeV (the scattered electron candidate) and the reconstruction in the central track- 
ing detector of the trajectories of exactly two charged particles (pion candidates) with opposite 
charges, transverse momenta larger than 0.1 GeV and polar angles confined within the interval 
20° < 9 < 160°. To reduce the background due to diffractive production of <p mesons, events 
with M KK < 1.04 GeV are discarded, where M KK is the invariant mass of the two hadron 
candidates when considered as kaons (no direct hadron identification is performed for this anal- 
ysis). In order to reduce both QED radiative corrections and background contributions in which 
there are unreconstructed particles, a cut E — p z > 52 GeV is applied. E—p z is the difference of 
the energies and the longitudinal momenta of the scattered electron (measured in the SPACAL) 
and the pion candidates (measured in the central tracking detector); it is expected to be close to 
twice the incident electron beam energy, i.e. 55 GeV, if no other particles have been produced 
except for the forward going system Y . To avoid backgrounds due to the diffractive production 
of systems decaying into two charged and additional neutral particles, all events are rejected in 
which a cluster, which is not associated with the electron or the two charged pion candidates, is 
reconstructed with polar angle larger than 20° and energy larger than 400 MeV (300 MeV) in 
the LAr (SPACAL) calorimeter. 

The p three-momentum is computed as the sum of the two charged pion candidate momenta. 
The variable Q 2 is reconstructed using the double angle method [|T0|]: 



2 _ AEj sinfl p (l + cosfl e ) 
V sm9 e + sm9 p - sin (9 e + 9 p ) ' W 

where Eq is the energy of the incoming electron and 9 e and 9 P are the scattered electron and 
p meson polar angles, respectively. The variable W is calculated using the Jacquet-Blondel 
method [[ft}: 

W 2 = ys-^- , with y = ^-^, (10) 
l-y 2E 

s being the square of the ep centre of mass energy and E p , p ZjP and p t>p being the energy, the 
longitudinal and the transverse momentum of the p meson, respectively. The electron transverse 
momentum is computed as 

2E -E p +p Z)P 

Pt,e — 7 !Ci /q\ • V iJ J 

The variable t' is then determined from the scattered electron and p momentum components 
transverse to the beam direction as 

t' ^ (Pt,miss) 2 = (Pt,e + Pt,p) 2 • (12) 



2 In the HI convention, the z axis is defined by the colliding beams, the forward direction being that of the 
outgoing proton beam (z > 0) and the backward direction that of the electron beam (z < 0). 
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The selected events are classified in two categories, corresponding to the absence or pres- 
ence of activity in the forward part of the HI detector. An event is classified in the "notag" 
sample when no signal above noise is detected in the PRT and the FMD, and no track and no 
LAr cluster with energy larger than 400 MeV is reconstructed with polar angle 9 < 20°. Con- 
versely, events are classified in the "tag" sample if a signal is observed in either the PRT or the 
FMD, or if a track or a LAr cluster with energy larger than 400 MeV is reconstructed in the 
forward part of the HI detector (9 < 20°). Q For the tagged events, a pseudorapidity interval of 
at least 2.2 unitsQ is required between the most backward track or LAr cluster with 9 < 20° (or 
the forward edge of the LAr calorimeter in the absence of track or cluster with 9 < 20°), and the 
most forward charged pion candidate. To first approximation, the notag and tag samples could 
be attributed to the elastic and the proton dissociative processes, respectively. However, elastic 
events fall in the tag sample when \t\ is large enough for the scattered proton to hit the beam 
pipe walls or adjacent material, leading to secondary particles which give a signal in the forward 
detectors. This effect becomes significant for t' > 0.75 GeV 2 . Conversely, proton dissociative 
events are classified in the notag sample for small masses, My < 1.6 GeV, or in the case of 
inefficiencies of the forward detectors. 

The uncorrected n + n~ mass distributions are shown over the extended mass region 0.3 < 
Mtttt < 1.3 GeV in Fig. [j], separately for the tag and the notag samples, for t' < 0.5 GeV 2 and 
for 0.5 < t' < 3 GeV 2 . Clear p meson signals are visible in all distributions. 



Monte Carlo simulations based on the DIFFVM program [ |T2| ] including QED radiation Ql3| ] 
are used to describe the production and decay of p vector mesons in elastic and proton disso- 
ciative scattering, and to correct the data for acceptance, smearing and radiative effects. The 
simulations include the angular distributions corresponding to the measurements of the present 
analysis for the matrix element (cos 9* distribution) and the Tq + 2r^ and rl + 2r| x com- 
binations ($ distribution). Other angular distributions and correlations are taken in the SCHC 
approximation, and the cos 5 parameter, which describes the interference between the longi- 
tudinal and transverse amplitudes, is taken from the elastic scattering measurement [JTp in the 
relevant Q 2 range. The exponential slope of the t distribution is b d = 7 GeV~ 2 for elastic 
scattering []lj] and b pd = 1.7 GeV -2 for proton dissociative scattering, values which describe 
well the t' distribution of the present data. The t'-integrated cross section ratio for proton disso- 
ciative to elastic scattering is taken as 0.75 in the present Q 2 range [|14|]. For proton dissociative 
scattering, the My spectrum is parameterised as da/dM Y oc 1/M y 15 (see JT5j]) and corrections 
are applied for the loss of events with large M Y values when particles of the dissociation system 
are reconstructed in the detector with polar angles 9 > 20° . All these parameters have been 
varied in the simulation as a part of the systematic error analysis. 

DIFFVM simulations have also been used for cu, <\> and p' background studies (see next 
section). In all cases, the t slopes are chosen to be b e i = 6 GeV~ 2 for elastic scattering and b p d = 
2.5 GeV~ 2 for proton dissociative scattering. The ratios of the proton dissociative to elastic 
channels, integrated over t', are 0.75. In the absence of measurements in electroproduction, the 
angular distributions for to, (except for — > K + K~ [Bp) and p' are treated as isotropic. 



3 In the case of proton dissociative scattering, this corresponds to an excitation mass of the target My 55 25 GeV. 

4 The pseudorapidity -q of an object with polar angle 9 is defined as 77 = — log tan(#/2). 
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Figure 1: Uncorrected it + it~ mass distributions for the selected events with 0.3 < 
M nn < 1.3 GeV. The shaded areas describe the uj and <p backgrounds and the hatched ar- 
eas the p' background, a) and c) show tag events with t'< 0.5 GeV 2 and 0.5 < t'< 3.0 GeV 2 , 
respectively; b) and d) show notag events. 
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2.2 uj, cf) and p' backgrounds 



Diffractive electroproduction of uo and mesons can fake p production through the decay chan- 
nels 



UJ — > 7T + 7T 7T° 



^ ^+^-^0 ; ^ ^o^o ^ (13) 

if the decay photons of the n° or the K\ meson are not detected. This happens if the deposited 
energy is associated with the charged pion tracks or does not pass the detection threshold in the 
detector. The p t imbalance of the event due to the loss of particles can then be interpreted as p 
production at large t', following eq. (|T2|). These background contributions, which are concen- 
trated below the selected p mass range (|7|), are estimated using the Monte-Carlo simulations. 
The ratios of the production cross sections uj I p and 4> I p are, for the present Q 2 range, taken 
as 0.09 and 0.20 [|], respectively. 



Another background reaction, particularly important for large t' in the selected mass range (|7|) 
is the electroproduction of p' mesonsf] decaying into two charged pions and two 7r°: 

fl -> p+TT^TT , p+ -> 7T+7T (+ C.C.) . (14) 

Again, the non-detection of the two 7r° mesons induces p t imbalance which fakes p production 
at large t'. No measurements exist in the relevant Q 2 range of the p' — ► 7r + 7r~7r°7r to p cross 



section ratio. The p' contribution is thus determined in section [O] from the data themselves, 
using the events with 0.5 < t' < 3 GeV 2 . 

It is important to recognise that the presence of backgrounds at large t' values affects dif- 
ferently the tag and the notag samples defined in section [ET} As mentioned there, genuine 
production of p mesons at large t', either due to proton dissociative or elastic scattering, usually 
gives a signal in the forward detectors and contributes mainly to the tag sample. In contrast, u, <p 
and p' background events, produced mainly at low t' but faking high t' p production, contribute 
to either the tag or the notag sample, depending on whether the proton dissociates and on the 
detector response. The ratio of the p signal to background at high t' is thus significantly higher 
in the tag sample than in the notag sample. 



2.3 Determination of the p' background 

In order to determine the p' background, a new variable, £, is introduced: 

^ Pt,miss ' Pt,p (15) 
\Pt ,miss 

For p' events produced at low t' and faking high t' p production, pt, m iss is due to the two missing 
7T° mesons and, in the present Q 2 range, is generally aligned along the n~) direction. This 

5 The detailed structure jl7| ] of the states described in the past as the p'(1600) meson is not relevant for the 
present study. The name p' is used to imply all vector meson states with mass in the range 1300-1700 MeV. In the 
simulations, the p' mass and width are taken as 1450 MeV and 300 MeV, respectively. 
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gives for the p' background a £ distribution peaking around +1, as shown in Figs. f|a-b; the 
same effect is found for to and production with the decay channels <JT3|) . In contrast, for 
genuine high t' p production, p t: miss is the transverse momentum of the scattered proton or 
baryonic system, leading to a flatter ( distribution, with maxima at —1 and +1. The ( and the 
$ distributions are strongly correlated: positive ( values correspond to $ angles close to 0° and 
360°, whereas for negative £, central $ values are selected. This is visible in Figs. 0c-d, which 
compare background events (for which £ is mostly positive) to p events (for which negative ( 
values dominate). 
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TAG EVENTS NOTAG EVENTS 
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Figure 2: Distributions of a-b) the ( and c-d) the $ variables of the selected events in the mass 
range (0) with 0.5 < t' < 3 GeV 2 , in the tag (a and c) and in the notag (b and d) sample. The 
shaded areas describe the uj and background as obtained from the simulations. The hatched 
areas correspond to the p' background and the open areas to the p contribution, as determined 
using the iterative fitting procedure described in section |23. 



An iterative fitting procedure is performed to estimate the p' background, whilst dealing 
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with this ( — $ correlation. After subtraction of the u and <\> backgrounds using the Monte 
Carlo simulations, the selected events in the mass range (0) are divided into four subsamples: 
tag events and notag events, separately with ( < and with ( > 0. Each of these four samples 
contains two contributions, due to genuine p production and to p' background. These eight 
contributions (eight unknowns) are determined through an overconstrained fit, performed using 
the MINUIT package []T8], to the numbers of events in the four samples (four measurements), 
under the following six constraints, obtained using the Monte Carlo simulations: 



i. two constraints describe the asymmetry of the ( distribution of p events (i.e. the ratio of 
the numbers of events with ( < and with £ > 0), separately for the tag and for the notag 
sample (see the open areas in Figs. [2]a and b, respectively); 

ii. similarly, two constraints describe the asymmetry of the ( distribution for p' events, sep- 
arately for the tag and notag samples (hatched areas); 

iii. the last two constraints, defined separately for £ < and for C > 0, are the probabili- 
ties for any p' event (elastic or proton dissociative) to be tagged; the ratio of the proton 
dissociative to elastic p' production cross sections is taken as 0.75 in the simulation. 



An estimate of the p' background is thus obtained, and its $ distribution is computed using 
the simulation. Relation @ is then fitted to the $ distribution in the tag sample, fully corrected 
for background, acceptance, smearing and radiative effects, to extract values of the spin density 
matrix element combinations Tq + 2r^ and Tq + 2r\ x . These values are fed back into the p 
simulation, leading to a modification of the simulated $ and hence ( distributions, which pro- 
vides new values for the constraints describing the asymmetry of the latter (see item i. above). 
The fitting procedure is repeated, and the iterative process converges after a few steps to stable 
background estimates, independent of the starting values of the spin density matrix elements in 
the p Monte Carlo simulation. 

Fig. || presents the £ and $ distributions of the selected events with 0.5 < t' < 3 GeV 2 , 
separately for the tag and the notag samples. They are well described by the superposition 
of the uj and <p background, the p' background and the p contribution, as determined from the 
iterative fitting procedure. The dominant background is found to be from p' production and, as 
expected, the backgrounds are larger in the notag sample and affect mainly the ( > region. 

This procedure thus provides an estimate of the p' / p cross section ratio for 0.5 < t' < 
3 GeV 2 . This estimate is extended to i! < 0.5 GeV 2 , under the assumptions quoted in sec- 
tion |2]T] for the t slopes and for the proton dissociative to elastic cross section ratio. 

The background contributions are shown in Figs. [TJa-d. After background subtraction, rel- 



ativistic Breit-Wigner functions, with the Ross-Stodolsky skewing parameter Q19| ] left free, are 
fitted to the fully corrected data, yielding p mass and width values in excellent agreement with 
expectations [ |T7| ] and good x 2 values. The data are thus very well described by diffractive p 
production with contributions of additional to, <p and p' backgrounds. The elastic cross section 
at low t agrees with previous measurements [[[]]. 
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2.4 Systematic errors 



In addition to the effect of varying the number of bins, the systematic uncertainties affecting the 
measurements described in section |3| are grouped into three classes: 

• Uncertainties in the amount and shape of the backgrounds 

The amount of uj and backgrounds is varied by ± 50 %. The uncertainty in the amount 
of p' background for t' > 0.5 GeV 2 is estimated by varying, for the fit procedure, the ( 
separation between the samples (at —0.4, —0.2, 0.2 and 0.4 instead of 0). The fractions 
of background events with ( < (see Figs. 0a-b) are also multiplied by 2 and 0.5. For 
t' < 0.5 GeV 2 , the p'/p cross section ratio is changed by ± 50 %. The shape of the 
background in the $ distribution (section EO) is modified by keeping the total amount 



fixed, but changing the fraction assigned to the two extreme $ bins by ± 50 % (see 
Figs. 0c-d). For the cos#* distribution (section |3T2"| ), the shape of the background is 



varied from flat to the same distribution as that of the p signal. The following model 
uncertainties in the background simulations are also included: the t slopes are varied 
(b e i = 6±1 GeV -2 andbpd = 2.5±1 GeV~ 2 ), the proton dissociative to elastic production 
cross section ratio is changed from 0.75 to 0.5 and to 1.0, and the mass and width of the 
p' meson are varied: M p > = 1450 ± 150 MeV and T p > = 300 ± 150 MeV. 

• Uncertainties affecting p production 

For the simulation of p meson production, the t slopes (b e i — 7 ± 1 GeV -2 and b p d = 
1.71q 7 GeV -2 ) and the proton dissociative to elastic production cross section ratio (0.75 
± 0.25) are changed, and the cross section dependences on Q 2 and W are varied within 
limits of the measurements in [JTp . Furthermore, the My spectrum as implemented in 
DIFFVM is varied from 1/M 2 / 15 to 1/M^ 85 and to 1/M^ 45 . 

• Uncertainties in the detector response 

The energy threshold for the detection of LAr clusters not associated to tracks is varied 
between 300 MeV and 500 MeV; the efficiencies of the PRT and FMD are modified 
within experimental errors; the measurement of the polar angle of the scattered electron 
is changed by ±0.5 mrad and the uncertainties in the trigger and the tracker efficiencies 
are included. 

For the measurements in section |3], the dominant systematic error is due to the uncertainty 
in the shapes of the backgrounds. 



3 t' dependences of spin density matrix elements 

3.1 <f> distributions and determination of (t*q + 2r^) and (r^ + 27^) 

Measurements of the spin density matrix element combinations r\ Q + 2rf x and r\ Q + 2r \ x are 
obtained from fits of eq. (Q) to the $ distributions, fully corrected for the presence of back- 
grounds and for acceptance, smearing and QED radiative effects, in the five t' bins shown in 
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Fig. For 0.5 <t'< 3.0 GeV 2 , only the tag sample is used, in view of the much larger back- 
ground in the notag sample (compare Figs, [ljc-d). Given the small backgrounds (Figs. [j]a-b) and 
in order to improve the statistical precision of the measurement, the tag and the notag samples 
are merged for if < 0.5 GeV 2 . The measurements are given in Table[l| and presented in Figs. 0a 
and |]b, together with previous measurements for t' < 0.5 GeV 2 in similar W and Q 2 ranges 
(Q 2 > 2.5 GeV 2 from ref. |p and Q 2 > 3.0 GeV 2 from ref. §]). 



Element 


t' (GeV 2 ) 


(0 (GeV 2 ) 


Measurement 


r 5 00 + 2r\ x 


t' < 0.08 
0.08 < t' < 0.2 
0.2 < if < 0.5 

U.O < t < l.U 

1.0 < if < 3.0 


0.037 
0.132 
0.320 
U. /UU 
1.620 


0.064 ± 0.012 ± 0.040 
0.087 ± 0.014 ± 0.017 
0.201 ± 0.014 ± 0.037 

U.lyo ± U.U1 / ± U.UiZ 

0.290 ± 0.023 ± 0.049 


r l m + 2r 1 1 1 


t' < 0.08 
0.08 <tf< 0.2 
0.2 < if < 0.5 
0.5 < if < 1.0 
1.0 < if < 3.0 


0.037 
0.132 
0.320 
0.700 
1.620 


-0.006 ± 0.025 ± 0.020 
-0.022 ± 0.027 ± 0.034 
-0.119 ±0.028 ±0.052 
-0.134 ±0.034 ±0.066 
-0.176 ±0.046 ±0.076 


r 04 
'00 


t' < 0.08 
0.08 < t' < 0.2 
0.2 < if < 0.5 
0.5 < t' < 1.0 
1.0 < t' < 3.0 


0.037 
0.132 
0.320 
0.700 
1.620 


0.678 ±0.013 ±0.015 
0.683 ±0.014 ±0.014 
0.662 ± 0.015 ± 0.025 
0.665 ± 0.019 ± 0.026 
0.708 ± 0.023 ± 0.035 



Table 1: Measurement of the spin density matrix element combinations r^ + 2r^, rg + 2r\ x 
and t[]q in five bins of t' . The first errors are statistical, the second systematic. 

Significant helicity non-conservation is observed in Fig. |]a for the combination Tq + 2r\ x 
(SCHC would imply a zero value of the combination). The r^ matrix element is proportional 
(see relations (Q)) to the product of the dominant non-flip amplitude T 00 and the T i amplitude, 
expected to be the largest helicity flip amplitude (see relations (Q)). In contrast, the rf x matrix 
element has a contribution from the product of the non-dominant non-flip amplitude Tn and 
the non-dominant single flip amplitude T 10 , and a contribution dependent on the double flip 
amplitude. The strong t' dependence of the r^ + 2r\ x combination is thus attributed mainly to 
the predicted [§-|8p \ftf dependence of the ratio of the T i to the non-flip amplitudes. Note that 
the t' dependence of the Tq + 2r\ 1 combination is not exactly oc Vf, as expected for the single- 
flip to the non-flip amplitude ratio, since it also involves the effect of the single and double-flip 
amplitudes in the denominator N of relations (0). 

The values for r^ + 2r\ x are shown in Fig. |]b. They are significantly different from zero 
and negative, which implies violation of SCHC. This is the first observation of a significant 
non-zero value of the r\ Q + 2r\ x spin density matrix element combination. The Tq element, 
which gives a negative contribution (see relations (Q)), is proportional to the square of the single 
flip amplitude T i. The r\ x element is proportional to the product of the non-dominant non-flip 
amplitude T u and the double flip amplitude. The sign of the combination thus gives information 
on the relative strength of the T iT^ and TuT^ products of amplitudes. It confirms that the 
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Figure 3: Normalised cross sections for p electroproduction in five bins in t'. The superimposed 
curves show the results of fits of a) relation @ and b) relation @. The error bars represent the 
statistical errors. 
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Figure 4: Measurement of a) rg + 2r^ 1 , b) rl + 2r{ x , c) as a function of t', for the present 
analysis and for refs. pip (labeled "HI elastic") and [@] (labeled "ZEUS M Y < 4 GeV"). The 
inner error bars represent the statistical errors, the full error bars include the systematic errors 
added in quadrature. The full curves correspond to the predictions of the model [|]] with param- 
eters extracted from the fit described in section |j. The dashed lines in a) and b) correspond to 
SCHC. 
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T 01 amplitude is significantly larger than the double flip amplitude in the present kinematical 
domain. The t' dependence of the 7q + 2r\ x combination is expected to be linear, up to effects 
of the single and double-flip amplitudes in the denominator N of relations (Q). 



3.2 cos 6* distributions and determination of 

The Tqq spin density matrix element is extracted from the cos 9* distribution using relation (TJ). 
Fig. |]b presents the fully corrected cos 6* distributions for five bins in t' . The normalisations of 
the uj, (f) and p' backgrounds are estimated as described above. The shape of the background is 
determined from the data and found to be flat. This is done by comparing the cos 9* distributions 
for the events with £ < and £ > 0, which differ in the data whereas they are predicted to be 
similar by the p Monte Carlo simulation. The difference is attributed to background. 

The extracted values of Tqq are presented in Table [j] and in Fig. |]c, together with previous 
measurements at low t' values in the same Q 2 and W range. No significant variation of Tqq with 
t' is observed. This is expected from relation (Q) in view of the predicted t' independence of 
the ratio |Tn|/|T o|, with small corrections from the other amplitudes. This observation implies 
that the slopes of the exponentially falling t distributions for the transverse and longitudinal 
s-channel helicity conserving amplitudes, T 00 and T n , are very similar. 



4 QCD description of the measurements 

Perturbative QCD calculations for vector meson electroproduction assume the factorisation of 
the non-perturbative from the perturbative contributions to the amplitudes. Collinear factori- 
sation has been demonstrated for longitudinal photons pop, but logarithmic singularities are 
manifest for the transverse photon polarisation when the fraction z of the longitudinal momen- 
tum carried by the quark approaches or 1 . Non-perturbative effects were suggested to damp 
these singularities []2T|], in which case perturbative calculations become problematic. However, 
as noted in p2|], these contributions cannot be large,f] and in the models [|]-|]] pQCD is expected 
to be valid also for the transverse amplitudes. 

In the models all amplitudes are proportional, in the leading log(<5 2 ) approximation, 
to the gluon density in the proton (except for non-perturbative contributions in the double-flip 
amplitude).[] More specifically, in the model of Ivanov and Kirschner (jjj, the gluon distribution 
in the proton xG(x, Q 2 ) is probed at the hard scale Q 2 = z(l — z)Q 2 < Q 2 /4. Following a 
suggestion of Martin, Ryskin and Teubner [|Z2|], the scale dependence of the gluon density in 
the leading log Q 2 approximation is parameterised as G(x, Q 2 ) = G(x, Ql) [Q 2 /QlY 1 ', where 
the gluon anomalous dimension 7 is taken as Q 2 independent. This permits the singularities as 
z — > 0, 1 to be avoided and factorisation is effectively (but not necessarily exactly) restored for 
the transverse amplitudes. 

6 Large non-perturbative contributions to the transverse amplitudes would imply a t slope similar to the non- 
perturbative photoproduction case, br — 10 GeV~ 2 , significantly larger than for the longitudinal cross section 



(6l — 6 GeV 2 ). This is inconsistent with the very weak t dependence of the Tqq matrix element (section 3.2 ) 
7 For large mass vector mesons and/or large Q 2 , skewed parton distributions should be used. 



17 



The absolute values and the Q 2 and t dependences of the four independent ratios^ of the 
amplitudes T\ VM \ are predicted by the model [|8|] with two independent parameters: the effec- 
tive gluon anomalous dimension 7, and the effective mass m of the incoming qq pair. These 
two free parameters are obtained from a fit to the t' dependence of the 15 measurements of spin 
density matrix elements in the present analysis.^ The fit gives an excellent description of the 
data: x 2 /ndf = 0.41 when the full errors are used and x 2 /ndf = 1.71 for statistical errors 
only. The fitted values of the parameters are 7 = 0.60 ± 0.09 and m = 0.58 ± 0.04 GeV. The 
errors are the quadratic combination of the statistical and systematic errors, the latter obtained 
by repeating the fits with the data shifted by each systematic uncertainty in turn. The dominant 
error comes from the uncertainty in the background shape. The results of the fit are shown as 
solid lines in Figs. |]a-c. As can be observed, the 11 low t' measurements of ref. [0]], which 
correspond to the same Q 2 and W ranges as the present data, are also very well described; their 
inclusion in the fit does not change the quantitative results significantly. 



According to the parton distributions in PDFLIB Q23|], the extracted value of 7 corresponds 
to the Q 2 evolution of the gluon density for Q 2 ~ 5.0 GeV 2 which is much higher than the aver- 
age (Q 2 ) < (Q 2 /4) ~ 1.3 GeV 2 in the data. The 7 parameter was introduced in the model to 
restore factorisation which is otherwise broken by end-point effects (z — > 0, 1) in the transverse 
amplitudes. It does not need to be strictly interpreted as describing the evolution of the gluon 
distribution in the proton at the specified Q 2 On the other hand, the disagreement may 
suggest that the model does not apply in the full Q 2 range of the present data. The parameter 
m = 0.58 ± 0.04 GeV is slightly below the p meson mass but belongs to a mass range where the 
quark pair is highly likely to recombine into a p meson (cf the parton-hadron duality arguments 
in [El). 



5 Conclusions 

A measurement has been performed of p meson diffractive electroproduction in the range 2.5 < 
Q 2 < 60 GeV 2 , 40 < W < 120 GeV and < f < 3 GeV 2 . The spin density matrix 
element and the combinations Tq + 2r\ x and 7q + 2r\ x have been measured as functions of f. 

No significant t' dependence is observed for r^. A significant violation of s-channel helicity 
conservation, increasing with t', is observed in the Tq + 2r\ x combination. It is consistent with 
a V¥ dependence of the ratio of the amplitude T ox to the non-flip amplitudes, T ox being the 
amplitude for the transition from a transverse photon to a longitudinal p meson. The Tq + 2r\ x 
combination is different from zero and negative; this is the first observation of a significant 
non-zero value of this combination. The sign gives information on the relative strength of the 
T 0X Tq X and T XX T\_ X amplitude products. Together with the Tq + 2r\ x measurement, it confirms 
that the T ox amplitude is relatively large in the present kinematical domain, and significantly 
larger than the double flip amplitude. 

A fit of the pQCD model of Ivanov and Kirschner to the present 15 measurements of 
spin density matrix elements gives a good description of the t' dependence of the data. The 

8 Under the assumptions of natural parity exchange and of purely imaginary amplitudes. 
9 The amplitude 2i_i, which contains a non-perturbative part in [|[] and is expected to be very small (see 
relations (Bb), can be set to zero or included as a free parameter in the fit without affecting the results. 
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value 7 = 0.60 ± 0.09 is obtained for the effective parameter describing the Q 2 dependence 
of the gluon density and m = 0.58 ± 0.04 GeV is extracted for the average effective mass 
of the incoming qq pair. Thus the data are broadly compatible with a diffractive mechanism 
based on the exchange of two gluons, with non-conservation of s-channel helicity occuring 
only when the photon longitudinal momentum is shared asymmetrically between the quark and 
the antiquark [||-|J]. 
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